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Abstract 

We re-investigate the question whether a light higgsino-like neutralino is a viable 
Dark Matter candidate. To this end we compute the dominant one-loop corrections to 
the masses of the higgsino-like states in the minimal Supersymmetric Standard Model 
(MSSM), due to loops involving heavy quarks and their superpartners. We also calculate 
analogous corrections to the couplings of higgsino-like neutralinos to Z and Higgs bosons. 
In the region of parameter space leading to high higgsino purity of the lightest neutralino, 
these corrections can change the expected relic density by up to a factor of five in 
either direction. We conclude that for favorable choices of soft supersymmetry breaking 
parameters, a state with more than 99% higgsino purity could indeed form all cold Dark 
Matter in the Universe. In some cases these corrections can also increase the expected 
cross section for LSP scattering off spinless nuclei by up to two orders of magnitude, or 
reduce it to zero. 


1) Introduction 

It has been known for more than ten years that the lightest supersymmetric particle (LSP) can 
be a good candidate for the missing “Dark Matter” (DM) in the Universe fli- In models with 
conserved “R parity” the LSP is absolutely stable. Searches for exotic isotopes f| then imply 
that it must be electrically and color neutral. Within the particle content of the Minimal 
Supersymmetric Standard Model (MSSM) this leaves us with two kinds of candidates, the 
lightest sneutrino v and the lightest neutralino Xv A combination of “new physics” searches 
at the CERN e + e~ collider LEP and direct DM search experiments excludes the sneutrino as 
viable DM candidate |2|. This leaves us with the lightest neutralino. 

In general the lightest neutralino y? is a superposition of the f/(l)y gaugino B (“bino”), 
the neutral 577 (2) gaugino W :i , and the two higgsinos h\ and h ° (with Y~ hi = — Y~ hn = —1/2). 
A gaugino-dominated state (photino or bino) has the right relic density f|, |5| to provide 
the missing Dark Matter if mj /m| 0 — (200 GeV) 2 , which points towards the same range 

of superparticle masses favored by naturalness arguments; here Ir stands for SU(2) singlet 
sleptons, whose exchange dominates the annihilation of bino-like neutralinos since they have 
the largest hypercharge. Originally it was thought |4|] that a sufficiently pure higgsino-state 
would also give an interesting relic density, if it is lighter than the W boson so that XiXi — 1 ► 
W + W~ is kinematically suppressed. The reason is that the coupling of a pair of higgsino-like 
LSPs to a Z boson becomes very small if the gaugino masses are much larger than the higgsino 
mass parameter. However, Mizuta and Yamaguchi [|| later realized that the standard estimate 
IH for the LSP density is not reliable in this case. The reason is that there are actually three 
higgsino-like states, two neutralinos and one chargino. The mass splitting between these states 
becomes very small as the gaugino masses become large. In such a situation one has to include 
“co-annihilation” between the LSP and these only slightly heavier states |7J. Note that the 
Z\ 1 X 2 and Wxixt couplings are large if the LSP is dominantly a higgsino. Co-annihilation 
therefore greatly reduces the estimate for the relic density of a higgsino-like LSP || [|. 

So far the discussion was based essentially on tree-level results (although QCD corrections 
were taken into account in the leading logarithmic approximation by using “running” quark 
masses and couplings ||). More recently, complete one-loop electroweak radiative corrections 
to the masses of the neutralinos and charginos have become available || . Using these general 
results, Giudice and Pomarol very recently pointed out U that loop corrections can quite 
significantly change the mass splitting between the three higgsino-like states, the dominant 
contribution coming from loops involving heavy quarks (f, b) and their superpartners. This is 
of some relevance for superparticle searches at LEP, since the heavier higgsino-like states will 
be quite difficult to detect experimentally if their decays only deposit a few GeV of visible 
energy in the detector |] 

Here we point out that these radiative corrections can also change the estimate of the 
LSP relic density quite dramatically, since the co-annihilation rate depends exponentially 
on the mass splitting between the higgsino-like states. Radiative corrections also change 
the decomposition of the LSP, which alters its couplings to gauge and Higgs bosons. These 
couplings are further modified by explicit vertex corrections. We present a full calculation of 
these three-point function corrections due to Yukawa interactions. Their effect on the relic 
density is relatively modest, but for negative sign of the higgsino mass parameter they can 

*In this case one can still search for events where the heavier higgsinos are produced in association with a 
hard, isolated photon. This signal should be viable even in the limit of almost perfect mass degeneracy, 
but the cross section is considerably smaller than for the simple pair production process. 
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change the cross section for LSP scattering off spinless nuclei by two orders of magnitude. 

The remainder of this article is organzied as follows. In Sec. 2 we describe the formalism, 
including one-loop corrections to the masses and relevant couplings of higgsinos. In Sec. 3 we 
present numerical results for the LSP relic density and its detection rate in a 76 Ge detector. 
Our estimate of the relic density includes a careful treatment of s —channel poles 0, as well 
as “sub-threshold annihilation” [|7j into W and Higgs boson pairs. Our estimate of the LSP 
detection rate includes the full set of contributions discussed in ref. []TT[ . In Sec. 4 we summarize 
our results and present some conclusions. Finally, in the Appendix we list expressions for 
three-point functions in the kinematical configurations of interest to us. 


2) Formalism 

In this section we describe the calculation of the one-loop corrections to the mass of the 
higgsino-like states, as well as to the couplings of the LSP to Z and Higgs bosons. We focus 
here on corrections from Yukawa couplings, which give the potentially largest contributions 
to the mass splittings || and couplings of interest to us. 


2a) Corrections to the Masses 

The corrections to the mass splittings can be understood as corrections to the chargino and 
neutralino mass matrices. Including one-loop corrections to the higgsino masses, these matri¬ 
ces can be written as 


M± = 


M 2 \/2M w sin/3 \ 
\/2 M w cos/3 ii + Sc ) 


(la) 


/ M l 

0 

—Mz cos/3 sin6W 
\ Mz sin/3 singly 


0 

m 2 

Mz cos/3 cos 6w 
—Mz sin/3 cos 9w 


—M z cos/3 singly 
Mz cos/3 cosOw 

$33 

—/l — A34 


M z sin/3 sin^vv \ 
—Mz sin/3 cos^w 

—/i — $34 
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(lb) 


Here Mi and M 2 are the supersymmetry breaking masses of the U(l)y and S77(2) gauginos, 
respectively, /i is the higgsino mass parameter, and tan/3 = (H®) /(Lf°) is the ratio of vacuum 
expectation values of the two neutral Higgs fields of the MSSM Jl2], |T3| . 

Note that we are only interested in states that are (almost) pure higgsinos. We do therefore 
not include contributions that are suppressed by both loop factors and small higgsino-gaugino 
mixing angles. This is why we only included corrections to the higgsino entries of the mass 
matrices. Further, since we focus on corrections due to Yukawa couplings, we only need to 
consider quark-squark loops. Their contributions are given by ||: 


<5c = 


32tt 2^ e { hbh t s[n ( 2d b) m t B o(Q,t,h) -B 0 (Q,t,b 2 

+h b li t sin(2 9 f )m b B 0 (Q, b , 3,) - B 0 (Q, b, t 2 ) (2a) 

+/J (hi sin 2 9~ b + h 2 t cos 2 0^j Bi(Q, t , (q) + (h 2 b cos 2 9~ b + h 2 sin 2 9-^j Hi(Q, t, b 2 ) 

+ (h 2 sin 2 9 b + h 2 cos^ 9ij i?i(Q, b , t±) + (ht cos 2 9 b + h 2 sin - 6ij B\(Q, b, t 2 ) |; 
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( 2 b) 

( 2 c) 

( 2 d) 


Here R 0 and B\ are two-point functions, for which we use the conventions of refs. || . Their 
first argument is the external momentum scale Q, and the second and third arguments are a 
quark and squark mass, for which we wrote the symbol of the corresponding fields in order to 


avoid double subscripts. The squark masses are eigenvalues of the t and b mass matrices [|T4 
which we write in the basis (qy, (Jr) following the notation of ref. [|T5[j : 


Mj = 
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m 2 t + mf r + (| — 


tL , V2 3 — Ow) cos(2 (5)M 2 Z 
—m t (A t + pcot/d) 
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(3b) 


Note that SU(2) invariance implies that the soft breaking terms appearing in the (1,1) entries 
of the mass matrices ([3a|) and ( |3T5j ) are equal. The lighter eigenstates are defined as q± = 
qLCOsOg + q/? sin Q~. Finally, hb and h t in eqs.(||) are the Yukawa couplings of the b and t 
quarks: 

gm b gm t 

\f2My/cos(3' t \f2MwAnf3’’ 

where g is the £77(2) gauge coupling, and the quark masses are to be taken at scale Q. 

As written, the corrections 5c and £34 are divergent. The fact that the divergence is the 
same for these two quantities provides a nontrivial check of our calculation. This divergence 
has to be absorbed by renormalizing the higgsino mass parameter fi. We have used the DR 
renormalization scheme, with renormalization scale taken equal to the external momentum 
scale Q. For consistency, the tree-level parameter /j in eqs.(|l]) then has to be interpreted as 
running mass taken at the same scale Q. In principle one has to diagonalize the matrices 
of eqs.([I|) at different Q = m^o or Q = m-± in order to compute the physical (on-shell) 
neutralino and chargino masses. However, since the Q —dependence of the corrections is quite 
weak, for our purposes it is sufficient to compute the mass matrices at fixed Q = \/i\. 

For our later discussion it is convenient to have approximate analytical expressions for the 
masses of the higgsino-like states as well as for the LSP eigenvector. In the for us relevant 
limit Mi, M 2 |p| the mass of the lighter chargino is approximately 

M\y sin(2/3) 


m-± ~ \fi C \ 


1 - 


HcM 2 


+ o(m 2 


- 2 \ 


( 5 ) 


where pc = A 4 + 5c] see eq. (|Ta|). 

In the same limit the two lightest neutralinos are approximately equal to the symmetric 
and anti-symmetric combination of the two higgsino current eigenstates h® and h 2 - Including 
terms up to first order in small quantities, their eigenvectors are given by: Q 

N h% ~ ( e l 5) ’ e 2 S) > + e ' 3 ’ _ es ) ; ( 6a) 

Hn our convention the neutralino eigenvectors are real, and we keep the signs of the eigenvalues. 
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with /i N = n + <5 34 . The upper (lower) sign in eqs.([7]a,b) holds for the symmetric (anti- 
symetric) higgsino state. The second equality in eq.flTcl) is valid only if one assumes the usual 
“unification condition” 

Mi = 1 tan 2 9wM 2 ~ 0.5M 2 . (8) 
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The masses of the higgsino-like eigenstates are given by 
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(9) 


where we have kept the signs of the eigenvalues, and the second equality again assumes eq.(||). 

If the second term in eq.(^) is larger than the loop corrections given by the third term, which 
is generally the case for M 2 < 1 TeV, the LSP will be the symmetric (anti-symmetric) higgsino 
state if /i is negative (positive). For small and moderate values of tan (3 this distinction is quite 
important, since the anti-symmetric higgsino-like state has larger gaugino components, see 
eqs.([/]a,b). Moreover, the mass splitting between the higgsino-like states also depends on the 
sign of /i. Assuming for simplicity eq.(|8]) to hold, one has || 


m { 0 — m^o ~ 
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— O-I-I — 0. 


5M 2 


'44 




m^o 

Xi 


M 2 1 

- [4 - sign (fj) sin(2/3)] - - (5 33 + <5 44 ) + S c - <5 34 . 


(10a) 


(10b) 


Note that the one-loop corrections increase (decrease) the mass splittings if <5 33 +<5 44 is negative 
(positive). As already pointed out in ref. 0 , this correction can be quite significant. In contrast, 
we found that the last term in eq.( |10b|) , Sc — 5 34 , is negligible in all cases. 


2b) Corrections to the Couplings 

At tree level the (} L coupling is proportional to the tree-level contribution to e 3 ; this 

coupling largely determines the annihilation rate of higgsino-like LSPs. However, at one- 
loop level one has to include the explicit vertex correction diagrams of Figs. la,b as well 
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as the (off-diagonal) wave function renormalization diagram of Fig. lc. We compute these 
corrections in the limit where the lighter neutralinos are exact higgsino; at the end we include 
mixing by simply multiplying the correction with the relevant higgsino component of the 
eigenstates. This procedure greatly simplifies the calculation. In this limit the only non¬ 
vanishing Z —neutralino coupling is the off-diagonal Z\iX 2 coupling, so we do not need to 
include any diagonal wave function renormalization diagrams. Note also that there is no 
vertex counter-term, since Yukawa couplings do not renormalize gauge couplings at one-loop 
level. Our procedure will give reliable results as long as the gaugino components e \, e 2 in 
eqs.(|j) are small. When these components become sizable or even dominant, our procedure 
may no longer give an accurate estimate of the loop corrections; however, in this case the 
loop corrections are in any case much smaller than the tree-level contributions to the mass 
splittings and couplings of interest to us, so that we again only make a small error. 



Fig. 1: Quark-squark loop corrections to the coupling of a pair of LSPs to a Z or Higgs boson. 
The LSP momenta k\ and k 2 point towards the vertex. Note that both senses of the “Dirac arrow” 
(flow of fermion number) have to be added, since the LSP is a Majorana fermion. There is also a 
diagram of type c) with a quark-squark bubble on the other neutralino line. There are two squark 
mass eigenstate with a given flavor. 


The diagrams of Fig. 1 can be described by the following effective ZxiXi vertex: 

3 9 


01 Zxx - 


87T 2 COS0W ■ 


(iv;^? + Nixaf) 7*75 + (A WW + A w®) (fcf + K) 75 


(ii) 

where N \3 and N \4 are the third and fourth component of the LSP eigenvector, and k\ and k 2 
are the momenta of the two neutralinos. We use the tensor decomposition of the three-point 
function as given in ref. [ 0 ; this form is convenient for the case of two equal external masses. 
The coefficients 5 a)P of eq.([TT|) can then be written as: 

= i c a,q + c v,q cos( 20 ,)] (m|o - k x ■ k 2 j c£(qi) + (m \0 + h ■ k 2 ^j C? (qi) 

+^*2(91) (0^0(91) — 2 C' 1 f ( 9 i)') + -mgCo(g 2 ) 
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+ C a < 
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-2I q sm 2 (29 q )C 2 {ch, q 2 )} 
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Here we have used the shorthand notation Ck(<ji ) = Ck(s,m^o,rn q ,m q ,m qi ), Ck(qi,qj) = 


Ck(s,m^o,rn qi ,rn qj ,m q ) and Bk(qi ) = Bk(mZ 0 ,m q ,m qi ). Recall that we use the B\ function of 


Xi 


refs. |Q], which differs from that of ref. [0 by an overall sign. Finally, the couplings in eqs.(|l2D 
are given by 


c = -I 3 ■ 
°a, q 2 3 ’Q j 


"■v, q g ± 6 q sin 9\y , 


(13) 


where I 3>q = ±1/2 and e q are the weak isospin and charge of quark 9, respectively. 

We note that one obtains a finite result only after summing over all three classes of diagrams 
and both squark eigenstates. Similarly, decoupling of degenerate heavy squarks (m qi ~ m q2 —> 
00) only holds after summation over all three diagrams and both squark eigenstates. However, 
the very last contribution to eq. ( |12a|) , is by itself finite and shows the proper decoupling 
behaviour. In fact, it closely resembles the corrections 5 33 and 5 44 of eqs.(|2]c,d). Indeed, these 
two terms come from the same two-point function diagrams, see Fig. lc. We can therefore 
include these terms either explicitly in or via the mass matrix corrections £33, <5 44 , where 
they change the quantity 63 given in eq.(|7c|); recall that this quantity determines the “tree- 
level” ZxiXi vertex fl2| : 


-pjU,tree _ • 9 

“*2cos 9 W 


7^75 


~ ±i- 


cos 9 


-— V2e 3 Y 


(N 2 l3 - N- 
75, 
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(14) 


where the (upper) lower sign is for the symmetric (anti-symmetric) higgsino-state, i.e. for 
negative (positive) ju. Together with eq.(|7c]) this sign ensures that the sign of the correction 
to the Zj(iXi coupling is independent of the sign of /i. This can also be seen from the 
last contribution to of course, keeping in mind that the LSP mass is always positive, 
independent of the sign of /i, see eq.(jP). We find that this term usually gives the dominant 
contribution to 5^; moreover, the 5^ are usually quite small. One can therefore get a 
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rough estimate of the size of the loop contributions from the diagrams of Figs. 1 by simply 
diagonalizing the mass matrix ([T5D, including the corrections £33 and £44, and using the “tree- 
level” vertex of eq. flUD . Of course, one must not include this correction both in the mass 
matrix and in 8^'. 

We also computed one-loop corrections to the LSP coupling to the neutral Higgs bosons of 
the MSSM. These couplings are not so important for the estimate of the relic density, unless 
2ro^o happens to be very close to or slightly lower than the mass of one of these Higgs bosons. 
However, the exchange of the neutral scalar Higgs bosons often gives the dominant contribution 
0.0 to elastic LSP-nucleus scattering. One therefore has to know these couplings quite 
accurately in order to make a reliable estimate of the event rate in various experiments that 
search for relic neutralinos. This is true both for direct detection experiments, which search 
for the recoil of nuclei struck by ambient LSPs, and for indirect detection experiments that 
search for neutrinos produced by LSP annihilation in the center of the Earth or Sun [|T^, [2J; 
the detection rate of the indirect search experiments is proportional to the rate with which 
ambient neutralinos are captured by the Earth or Sun, which in turn is proportional to the 
LSP scattering cross section off ordinary matter. 

At tree-level a pure higgsino state has no couplings to Higgs bosons; these couplings 
originate from the Higgs-higgsino-gaugino interactions in the supersymmetric Lagrangian. 
This also implies that the wave function renormalization diagram of Fig. lc does not contribute 
liere.f] We therefore only have to evaluate the explicit vertex corrections of Figs. la,b. 

Their contribution can be described by the effective vertices 
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where (j) stands for the light neutral scalar h° or the heavy neutral scalar H° and A is the 
pseudoscalar Higgs boson. The coefficients in eqs. (p~5|) can be written as 
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1 There are wave function renormalization diagrams where an external higgsino is converted into a gaug- 
ino, which then couples to the Higgs boson and the second higgsino. This gives contributions of order 

j^g 2 htrrit/M 2 , which can be interpreted as 0(h 2 ) corrections to the tree-level coupling, which is of or¬ 
der gM\y/M 2 . However, unlike our corrections ^33 and 644 this only corrects an entry in the neutralino mass 
matrix that is already nonzero at tree-level, and will therefore change the LSP couplings to Higgs bosons 
by at most a few percent. In our case these corrections are further suppressed since we are interested in a 
higgsino-like LSP, which (for tan/3 7 ^ 1) implies Mi, M 2 3> diagrams with internal gaugino lines are then 
suppressed by the large gaugino masses. 
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Here we have used the same notation for the arguments of the C functions as in eqs. fll^) . The 
coefficients and describe the Higgs couplings to quarks; they are given by [HI : 
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where a is the mixing angle of the neutral scalar Higgs bosons jnj. Finally, the coefficients 
M and c ~\2 describe the couplings of one Higgs boson to a pair of squarks; they are given 

by m 
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The A q also appear in the squark mass matrices of eqs.(^), and I 3:q and e q again refer to 
the weak isospin and charge of quark q. Finally, the Yukawa couplings h q have been defined 
in eq.([I|). 

Note that the pseudoscalar Higgs boson has no couplings to two equal squark eigenstates 
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We note that in this case the diagrams of the type shown in Fig. la are finite by themselves 
once one has summed over both squark eigenstates, and the diagrams of Fig. lb are separately 
finite for each combination of squarks in the loop. Notice also that this last class of diagrams 
is proportional to the Higss-squark-squark couplings, which receive contributions from the A q 
parameters; these couplings can become very large 0. 


As noted earlier, we use the C functions of ref. |H|. However, there is a technical complica¬ 


tion. When estimating the LSP relic density we need to evaluate these functions at s = 4m~ 0 , 
whereas LSP-nucleus scattering cross sections probe these functions at s ~ 0. The expressions 
for the higher C functions given in Appendix C of ref. ]16[ contain apparent divergencies in both 


these limits. We stress that the loop functions themselves remain well behaved as s 


4 m 2 


x? or 

s —> 0; the apparent divergencies in the expressions of ref. |16| therefore all cancel. In fact, even 
the standard expression for the scalar three-point function Cq contains apparent divergencies 
in the kinematical configurations of interest to us. In the case of C 0 the necessary cancella¬ 
tion between different terms can still be accomplished numerically by slightly increasing or 
reducing s. However, for the higher C functions these cancellations become quite delicate. 
We therefore re-evaluated the relevant Feynman integrals for the two cases of interest to us. 
In both limiting situations the C functions can be expressed as combinations of two-point 
(B) functions; all coefficients are now finite. The relevant expressions are collected in the 
Appendix. 


3) Results 


We are now in a position to present some numerical results. We focus on light higgsino-like 
states, m^o < M\y , since heavier LSPs have very large annihilation cross sections into W 
and Z pairs 01- A heavy higgsino therefore only makes a good cold Dark Matter (CDM) 
candidate if its mass exceeds 0.5 TeV. This is already uncomfortably heavy for “weak scale” 
supersymmetry; for example, assuming gaugino mass unification, the gluino mass has to be 
larger than 3 TeV in such scenarios. In fact, the annihilation cross section into W + W~ final 
states is so large that it can be relevant even if the LSP mass is a little below M]y. Such “sub¬ 
threshold annihilation” can occur since at freeze-out the LSPs still have significant thermal 
energy. We include this effect for W + W~ and also h°h° pairs in our estimate of the LSP relic 
density, using the formalism developed in ref.0. We also use a careful treatment of s —channel 
poles (Z and Higgs exchange diagrams); as pointed out in ref.Q, the standard expansion in 
the LSP velocity [II] breaks down in the vicinity of such poles. We use the numerical method 
developed in ref. [[2T |. 

In order to illustrate the effects of the loop corrections on the LSP couplings to Higgs 
bosons, we also present results for the LSP counting rate in an isotopically pure ,6 Ge detector, 
assuming a fixed local LSP mass density of 0.3 GeV/cm 3 |32[ and a velocity dispersion of 
320 km/sec. Nuclear effects are described by a Gaussian form factor, with a nuclear radius 
of 4.1 fm [Q. Of course, we could just as well have used any other spinless isotope. The 
scattering rate due to spin-dependent interactions is affected by the loop corrections to the 
ZxiXi coupling, but this correction is usually somewhat smaller than that to the LSP-Higgs 
couplings. Note that the total scattering cross-section off heavy nuclei is usually dominated 
by the spin-independent contribution even if the nucleus in question does have non-vanishing 

M. 


spin 


Since we are interested in scenarios with rather light LSP, we have to be careful not to 
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violate any experimental bounds. The most relevant constraints on the parameters appearing 
in the tree-level neutra lin o mass matrix conies from chargino searches at LEP |25|]. Unfortu¬ 
nately these bounds are not entirely straightforward to interpret in our case, since the standard 
set of experimental cuts used to suppress SM backgrounds becomes quite inefficient in sce¬ 


narios with small Am,- = m~± 

* Xi 


m^o. Note also that the cross section for the production of 


higgsino-like charginos is smaller than for gaugino-like states. We interpret the LEP bounds 
as requiring 

' 75 GeV, 


. . . „ „„ . , A rriy > 10 GeV 

Xf - 1 45 GeV, Am, < 10 GeV 


( 20 ) 


The second bound comes from the measurement of the total width of the Z boson m- and 
thus holds for any value of A m^. We are aware that this parametrization of the LEP search 
constraints is only a crude approximation, but it should be sufficient to illustrate the effects 
of the loop corrections. 

We have seen in the previous section that these corrections depend on the details of the stop 
and sbottom mass matrices. In particular, the corrections £33 and 8 . 44 to the neutralino mass 
matrix are proportional to sin(20q), see eqs.(|2]c,d); these corrections also vanish in the limit of 
equal masses for squarks of a given flavor. The combination of these two properties means that 
the corrections depend sensitively on the size of the off-diagonal entries of the squark mass 
matrices (|3j). Moreover, the potentially largest correction to the LSP-Higgs coupling, coming 
from the diagram of Fig. lb, directly depends on the A—parameters appearing in the squark 
mass matrices. Third generation squarks also contribute to other loop processes. This imposes 
some constraints even on combinations of parameters where all squark mass eigenstates lie 
well above the direct experimental search limits |26|, 27 . 

In ref. 0 the t — b loop contribution to the electroweak p parameter was emphasized. 
However, given that a “new physics” contribution ~ 3 • 1CT 3 is not excluded by present 


data ]27], we find that other loop corrections lead to stronger constraints. I 11 particular, loop 
corrections to the mass of the light neutral Higgs scalar h° turn negative when A t becomes 
too large |^S|. One important constraint therefore comes from searches for the MSSM Higgs 


bosons at LEP 27 


The constraints we have discussed so far do not depend on the masses of the other squarks, 
and are therefore quite model-independent. If we make the additional simplifying assumption 
that all explicitly supersymmetry breaking diagonal squark masses (m^, m^ R and 'm~ bR in 
eqs.(||), and analogous quantities for the first and second generation squarks) are equal at the 
weak scale, we find that the strongest constraint on the parameters of the stop mass matrix 
comes from the recent CLEO measurement [^UJ of the branching ratio for inclusive b —> sy 
decays: 

1 • 1CT 4 < B(b —> sy) < 4 • 1CT 4 . (21) 

Since we are studying scenarios with rather light charginos, chargino-stop contributions to this 
partial width can be quite large |30|; they can be of either sign, depending on the signs of p 
and A t . However, the resulting constraint is more model-dependent: If one allows some non¬ 
universality of soft-breaking squark masses, one also gets contributions from gluino-squark 
and neutralino-squark loops f3(|, the size of which depends strongly on the details of the 
entire three generation squark mass matrices. For defineteness we will stick to a scenario with 
exactly universal soft breaking squark masses, and with A t = Ab = A, with the understanding 
that the constraints that result from imposing the bounds ( |2T| ) can be relaxed in slightly more 
general models without significantly changing the loop corrections to the masses and couplings 
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of higgsino- like states. 

In Fig. 2 we show the dependence of various quantities relevant to our subsequent analysis, 
normalized such that they can be plotted to a common scale. We fixed M 2 = 350 GeV and 
H = —70 GeV, which means that the LSP is a more than 99% pure higgsino; we define the 
higgsino fraction as 1 — gaugino fraction = 1 — (ef + e 2 ), see eqs.(||). We chose tan/3 = 1.5 
so that the top Yukawa coupling is close to its upper bound, if one requires it to remain 
perturbatively small all the way to the GUT scale; on the other hand, b — b loops are essentially 
negligible for such a small value of tan/3. We took a very large mass (1.5 TeV) for the 
pseudoscalar Higgs boson, since this maximizes m^o, and hence minimizes the impact of the 
LEP Higgs search bounds. This also means that charged Higgs boson loop contributions to 
the b —> sy partial width are negligible. Our choice of 430 GeV for the common soft breaking 
squark mass is again motivated by our desire to maximize the size of the loop effects, given 
the experimental constraints discussed above. Increasing trig for fixed Ajrriq would reduce the 
ratio of physical stop masses, which leads to reduced t — t loop corrections. On the other 
hand, we cannot increase A/rriq beyond the limits shown in Fig. 2 without violating some 
experimental bound. Finally, here and in the subsequent figures we assume gaugino mass 
unification, eq.(§§). 

The curves in Fig. 2 terminate at values of A where m/jO falls below the LEP bound of 
about 62 GeV; note that h° is essentially indistinguishable from the single Higgs boson of 
the SM if m 2 A M§. The two dotted curves show a “high” and “low” theoretical estimate 
for B(b —> sy), scaled up by a factor 10 4 . Our estimates are based on a leading order QCD 

m 


analysis 


which has substantial scale uncertainties 


the band in Fig. 2 corresponds 


to varying the renormalization scale between 2.5 and 10 GeV, and also includes uncertainties 
from CKM matrix elements etc.| Notice that the “low estimate” can be zero. This happens 
if the contribution from sparticle loops is larger than that from the standard t — W loops 
and has opposite sign, reversing the sign of the complete matrix element at scale My y or m t . 
Renormalization group effects give another contribution from tree-level W exchange due to 
operator mixing; this contribution is not sensitive to any “new physics”. In the SM this term 
has the same sign as the loop matrix element at scale Myp, leading to a large QCD enhance¬ 
ment factor, but in the MSSM these two contributions can cancel. For fixed renormalization 
scale this cancellation only happens at specific points of SUSY parameter space, but perfect 
cancellation becomes possible for an entire range of parameters if the renormalization scale is 
allowed to vary. 

We note that for the given sign of /i, the branching ratio for inclusive b —► sy decays tends 
to be below (above) the SM prediction if A is negative (positive). In order to be conservative, 
we only exclude combinations of parameters where the “high” theoretical estimate is below the 
lower bound, or the “low” estimate is above the upper bound, given in (|2l|) . For the parameters 
of Fig. 2 this translates into the constraint —2.7 < A/rriq < 2.65, which is only slightly stronger 
than that resulting from Higgs searches at LEP. Within this region, 5pf b < 2.2 • 10~ 3 . 

The solid curve in Fig. 2 shows Ar/y (divided by 5 GeV to fit the scale). The tree-level 
prediction for this quantity for the given choice of M 2 , p and tan/3, 14.5 GeV, is very close 
to the loop-corrected value for A — 0. We see that the corrections can either increase or 
decrease the chargino-LSP mass splitting by about 4 GeV before one gets into conflict with 
the constraint (|2l|) . In this case the loop corrections therefore only amount to at most 30%; 

sy) in the framework of the 


Very recently an almost-complete next-to-leading order calculation oiB(b 
MSSM has appeared p2|; their result falls within our band. 
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however, we will see below that this suffices to change the prediction for the LSP relic density 
quite dramatically. 

Finally, the long and short dashed curves in Fig. 2 show the (re-scaled) couplings of 
the LSP to Z and h° bosons, respectively; gzx°x° defined as the axial vector coupling at 
s = 4m| 0 , while g h ^ o^o is defined at s = 0. As in case of Am*, the tree-level predictions for 
these quantities are very close to the loop-corrected values at A = 0. We see that the relative 
variation in LSP —Z coupling is larger than that in Am^, when A is varied over its allowed 
range. Note also the positive correlation between these two quantities, which reinforces the 
correlation between small Aand small gzx°x° that holds for higgsino-like LSPs at tree- 
level, see eqs.(|7cD and (|10b| ) . A similar correlation also holds for the loop-corrected coupling 
of the LSP to the light scalar Higgs boson, as shown by the short dashed curve. However, in 
this case the tree-level prediction is very small, g h ^ .o^o ~ 7.6 • 1CT 3 . This can be understood 
from eqs. (0a, b) and the general expression for this coupling given in ref. [|li|] : 


3h°x ?Xi;tree — 


\/2 


( (S,A) , ( S,A)\ , . , >. 

ge 2 — g e 1 I (sina ± cosa) 


- |(sin/?^cos/3) 2 ^|^, 

5 iKZ 2 


( 22 ) 


where g' = gtand w is the U( l)y gauge coupling, and the upper (lower) signs again hold for 
the symmetric (anti-symmetric) higgsino; in the second step we have used the “unification 
condition” (||) as well as the relation a ~ (5— |, which holds for m 2 A M§. For the quite small 
value of tan/? used in Fig. 2 this gives a strong cancellation in the coupling of the symmetric 
higgsino-like state, which is the LSP for fi < 0. As a result, the one-loop correction can easily 


dominate over the tree-level contribution (^2|). This leads to the behaviour shown in Fig. 2, 
where the coupling changes sign at A ~ —0.8 rriq. 

In Figs. 3a-c we show the chargino-LSP mass splitting, the LSP relic density, and the LSP 
detection rate in a 76 Ge detector as a function of the gaugino fraction of the LSP eigenstate. 
We have again chosen rriq = 430 GeV, tan/? = 1.5 and a large value of rriA■ Note that, unlike 
in Fig. 2, the physical LSP mass has been kept fixed in Figs. 3; the value of 70 GeV chosen 
here is close to that which maximizes the estimate of the relic density. Since rri^o is kept fixed, 
both M 2 and the tree-level parameter g(\g\) vary along the curves; e.g., M 2 lies between about 
150 GeV and 1 TeV, with larger values of M 2 corresponding to smaller gaugino fractions, see 
eqs. 0a, b). In order to maximize the loop effects we have also varied the A—parameter slightly. 
In the region of relatively small M 2 , i.e. large gaugino fraction, the light chargino is somewhat 
heavier, as shown in Fig. 3a; this reduces the absolute size of the t — loop contributions to 
the b —> sy decay amplitude for fixed A, which in turn allows us to go to slightly larger values 


of |A| without violating the bounds (|2l|). 

We show three curves in each of Figs. 3. The dotted curves labelled “no loops” have been 
obtained by switching off the loop corrections discussed in Sec. 2. However, we keep quark- 
squark loop contributions to the mass matrix of the scalar Higgs bosons fh| [28]|, as well as q 
loop contributions to the h°gg coupling and q — q loop contributions (box diagrams) to the 
LSP-gluon coupling [34, |HJ]. These corrections depend only weakly on the sign of A, however. 
On the other hand, the signs of the corrections discussed in Sec. 2 are essentially fixed by the 
sign of A, as shown in Fig. 2. In Figs. 3 we therefore show results both for positive (solid) 
and negative (dashed) A, keeping |A| fixed. Since we chose parameters close to those that 
maximize these loop corrections, the band between the solid and dashed curves in Figs. 3a,b 
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roughly indicates the range that can be covered by changing the parameters of the squark mass 
matrix, for fixed values of the parameters appearing in the tree-level chargino and neutralino 
mass matrices. 

The results of Fig. 3a show that loop corrections can change the chargino-LSP mass 
splitting by about three to four GeV in either direction, as already indicated in Fig. 2. Note 
that the absolute size of this correction is almost independent of the gaugino fraction. This can 
be understood from eq. (| 10 b|) , which shows that the tree-level and loop-induced contributions 
to Aare independent of each other as long as the LSP is a higgsino-like state. 

The results of Fig. 3b show that in the region of high higgsino purity the relatively mod¬ 
est loop corrections to A my can change the estimate of the LSP relic density by more than 
a factor of five. The reason is that here the relic density is essentially determined by co¬ 
annihilation processes |J, which depend exponentially on the mass splitting [[7|. Our calcu¬ 
lation includes Xixf co-annihilation into //' and IMy final states through W exchange, and 
X 1 X 2 co-annihilation into // final states through Z exchange, where / stands for any SM 
fermion other than the top quark. As is usually done, we show results for the LSP mass 
density in units of the critical or closure density, Qy = p^o /p c , multiplied with the square of h, 
the ffubble constant in units of 100 km/(sec-Mpc); a conservative range for h is 0.4 < h < 0.9, 
with recent measurements clustering around 0.5 — 0.6. One needs fl^h 2 > 0.02 — 0.03 if the 
LSP is to form the bulk of the galactic Dark Matter haloes, and Vt^h 2 > 0.15 if the LSP is to 
form all CDM in models |36|] with mixed hot and cold Dark Matter. The results of Fig. 3b 
show that, if A is large and positive, a 99.9% pure higgsino state can form galactic haloes, and 


a 99.5% pure higgsino might form all CDM. On the other hand, if A is large and negative, 
one will need at least 1% gaugino fraction, corresponding to 61,2 ~ 0.1, even for the LSP to 
be able to form galactic haloes. 

Note that the curves in Fig. 3b cross over in the region where the gaugino fraction exceeds 
several percent. Here the mass splitting between the higgsino-like states becomes so large that 
the relic density is again determined by the usual X 1 X 1 annihilation processes, which in our 
case mostly proceed through virtual Z exchange. We saw in Fig. 2 that the loop corrections 
increase (reduce) the ZxiXi coupling if A is positive (negative). As a result, the curve for 
A > 0 reaches its maximum already at a rather small gaugino fraction; in this case the LSP 
can form all CDM if 180 GeV < M 2 < 340 GeV. In contrast, the curve for A < 0 reaches its 
maximum at larger gaugino fraction; here the LSP can from all CDM only if M 2 falls in the 
narrow window between 160 and 195 GeV. 

In Fig. 3c we show estimates for the LSP detection rate in , 6 Ge, ignoring possible energy 
thresholds and assuming a fixed local LSP mass density. In this case the loop corrections 
discussed in Sec. 2 can increase the tree-level result by more than two orders of magnitude! 
This is largely due to the small tree-level value of f° r ^ ie given case of a symmetric 

higgsino-like state, see eq.(|22|). The turn-over in the region of sizable gaugino fraction is 
caused by mixing with the bino-like neutralino; note that near the end of the curves shown 
in Figs. 3, Mi and /i are already quite close to each other, so the expression flU) for ei is no 
longer reliable. 

We emphasize that in this case the band between the solid and dashed curves is not a 
good estimate of the variation of the expected counting rate when the parameters of the stop 
mass matrix are varied, since for large |A| the loop corrections dominate over the tree-level 
contribution to g h ^ 0 ^. 0 , as shown in Fig. 2. The total scattering rate can be made to vanish 
exactly for moderately negative values of A; this is related to the change of sign of 57 ^ 0^0 
observed in Fig. 2. Note that the scattering rate can vanish even for quite moderate values of 


13 








all sparticle masses. This illustrates that it is impossible to give a strict lower bound on the 
expected LSP detection rate even within the MSSM. Of course, there is no a priori reason for 
such a cancellation between tree-level and one-loop contributions to occur; indeed, over most 
of the parameter space the loop corrections increase the expected event rate. However, even 
the most optimistic estimate in Fig. 3c is still several orders of magnitude below the sensitivity 
of present experiments ][37|| . 

Recall that we do not rescale 


38] the event rate in regions of parameter space leading to 


a very small LSP relic density; had we done so, most of the curve for A < 0 would have been 
below the one for positive A, as can be seen from Fig. 3b. Finally, we saw in Fig. 2 that for fixed 
|A|, the absolute size of g h ^ o^o is somewhat smaller for A < 0 than for A > 0. We nevertheless 
find a slightly larger scattering rate for A < 0, partly because this also gives a slightly smaller 
value for m^o; the h° exchange contribution to the LSP-nucleon scattering matrix element 
scales like g h ^o^/m\ 0 - Destructive interference with various squark loop diagrams [|TT[ also 
plays a role here. 

In Figs. 4a-c we show results similar to those of Figs. 3, but for positive sign of the higgsino 
mass parameter p. The choices for the other parameters are very similar to those in Figs. 3, 
except that A is now fixed along each curve. For this rather small value of tan/3, flipping the 
sign of p has quite dramatic effects, as already anticipated in our discussion in Sec. 2. In 
particular the gaugino fraction of the LSP for fixed values of M 2 and |p| has become much 
larger. Conversely, one has to go to much higher M 2 in order to achieve a given level of 
higgsino purity; in Figs. 4, M 2 varies between about 0.3 and 1.3 TeV. This also implies that 
for given higgsino purity the chargino-LSP mass splitting is smaller for p > 0 than for p < 0, 
see eq. (|10b|) . On the other hand, Fig. 4a shows that the size of the loop contributions to this 
mass splitting is essentially independent of the sign of p, as long as the gaugino fraction is 
small. Note that for positive p, Acan be below 10 GeV for a gaugino fraction as large as 
12.5% (e li2 ~ 0.3); this will have ramifications for chargino searches at LEP 0. 

The smaller Afor fixed gaugino fraction also implies a greatly reduced relic density, due 
to enhanced co-annihilation rates. This is illustrated in Fig. 4b. The increased importance 
of co-annihilation also helps to explain why the curves in this figure do not cross, in contrast 
to those in Fig. 3b. Another reason is that in the region of sizable gaugino fraction, £3 is 
smaller for positive p; this is due to non-leading, (P(Mf 2 ) terms not included in eq.flTcj), 
which become quite important when the gaugino fraction exceeds several percent. Indeed, for 
the largest gaugino fractions shown in Figs. 4, our treatment of the loop corrections to the 
couplings of the LSP may no longer be entirely reliable, as discussed in Sec. 2 ; however, as 
anticipated in the same discussion, the relative importance of the loop effects decreases with 
increasing gaugino fraction. 

This is also true for the estimated LSP detection rate shown in Fig. 4c. For the smaller 
gaugino fractions shown in this figure, we find that the loop corrections can change the estimate 
by a factor of about two in either direction, whereas for large gaugino fraction the loop effects 
amount to at most 30%. We note again that, had we re-scaled the event rate for scenarios with 
small LSP relic density, the loop effects would have been even more important in the region of 
high higgsino purity. Finally, note that in the region where the gaugino fraction exceeds 5% 
the counting rate in Fig. 4c exceeds that in Fig. 3c by almost an order of magnitude. This is 
due to the much larger tree-level value of g h ^ o^o for the anti-symmetric higgsino-like state, see 
eq. (| 22 |) . Since now the tree-level value exceeds the loop corrections, we find smaller (larger) 
counting rates for negative (positive) values of A. 

As a final illustration of the effects of the loop corrections to the masses and couplings of 
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higgsino like states, we show in Figs. 5a-c the “geography” of the well-known (M 2 , /i) plane in 
the region M 2 \fjb\, fi < 0. In Fig. 5a these loops have been switched off, while in Figs. 5b,c 
they have been included with A = —2.5 rrig and A = +2.5mg, respectively. In each case the 
region to the right of the solid line is excluded by the LEP chargino search limit (^). The 
long and short dashed curves are contours of constant LSP relic density Q^h 2 = 0.025 and 
0.1, respectively. The remaining lines are contours of constant LSP detection rate in a ' 6 Ge 
detector, measured in events/(kg-day); as before, we have assumed a fixed local LSP density 
when calculatimg the counting rate. 

We see that, depending on the sign of A, loop corrections to the chargino and neutralino 
mass matrices can signhcantly reduce (Fig. 5b) or increase (5c) the size of the region that is 
excluded by chargino searches at LEP; this is a direct result of the change in Adepicted 
in Fig. 3a. Similarly, the loop corrections can increase or decrease the region where the LSP 
is a good CDM candidate; recall that A 0.025 is required if LSPs are to form the bulk 
of galactic Dark Matter haloes. Finally, we again observe a dramatic change of the expected 
LSP detection rate due to radiative corrections. Note that this rate changes only very slowly 
in Figs. 5b,c, whereas at tree-level one expects a significant dependence on M 2 , see Fig. 5a. 

We note in passing that results similar to those displayed in Figs. 3 and 5 can also be 
obtained in the framework of a recently proposed [[kj model with non-unihed gaugino masses 
and a higgsino-like LSP. This model attempts a supersymmetric interpretion of an event 
with an e + e _ pair, two hard photons, and missing transverse momentum px reported by 
the CDF collaboration fLl. The prospects for detecting relic neutralinos in this model have 


recently been studied in refs.[j4T|, using tree-level results for the neutralino mass matrix and 
LSP couplings. We expect that loop corrections of the type discussed here can modify these 
estimates significantly. However, searches for additional events with two hard photons and 

H: 


missing p T failed to find additional candidates 


casting doubt on any supersymmetric 


interpretation of the single anomalous event. We therefore do not study this model in any 
further detail. 


4.) Summary and Conclusions 

In this paper we have presented a calculation of loop corrections involving Yukawa couplings 
to the masses and couplings of the higgsino-like states of the MSSM. We have found these 
corrections to be very sensitive to the size and sign of the soft supersymmetry breaking A 
parameters. If |H| is large, the one-loop prediction for the difference of the chargino and LSP 
masses can differ by up to ~ ±4 GeV from the tree-level estimate; the loop correction to the 
X 2 — Xi mass difference is about twice as big. Combinations of parameters leading to even 
larger corrections lead to conflicts with the measured value of the branching ratio for inclusive 
b —*■ S 7 decays, and/or with the negative outcome of searches for Higgs bosons at LEP. We 
also found that for negative sign of the higgsino mass parameter //, one-loop corrections to the 
LSP coupling to Z and Higgs bosons can be comparable to or even larger than the tree-level 
contributions already for quite moderate gaugino masses, M 2 > 200 GeV. 

We have illustrated the importance of these loop corrections by computing their effect on 
the estimated LSP relic density and on the direct LSP detection rate, assuming the LSP to be 
higgsino-like. The relic density is in this case often determined by co-annihilation processes, 
the rate of which depends exponentially on the mass splittings between the higgsino-like states. 
Yukawa loop corrections can therefore change the tree-level prediction by a factor of ~ 5 in 
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either direction. This re-introduces a state with more than 99% higgsino purity as a viable 
cold Dark Matter candidate, if the corrections to the mass splittings are near the upper end of 
their allowed range. If /i < 0, the effect of loop corrections on the estimated LSP counting rate 
is even more dramatic: The predicted rate might increase by two orders of magnitude, but it 
might also be exactly zero (for spinless nuclei), even if all sparticle masses are in or below the 
few hundred GeV range. Clearly effects of this size have to be included in any quantitative 
analysis of the properties of higgsino-like Dark Matter. 

We conclude with some remarks regarding the viability of models with higgsino-like LSP. 
Within the framework of minimal supergravity models [|43|, which assume universal scalar 
masses as well as unified gaugino masses at the Grand Unification scale, a higgsino-like LSP 
is possible only if tan 2 /? 1, and if scalar soft breaking masses are significantly larger than 
gaugino masses. Since gaugino masses in turn must be considerably larger than \/i\ for the 
LSP to be higgsino-like, naturalness arguments favour a very light LSP in such a scenario. 
On the other hand, the parameter space leading to a higgsino-like LSP opens up considerably 
if one allows the sparticle spectrum at the GLIT scale to be non-universal. In particular, 
the predicted value of // can be reduced either by giving larger soft breaking masses to the 
Higgs bosons than to third generation squarks, or by reducing the gluino mass compared 
to the masses of the electroweak gauginos. We therefore conclude that a higgsino-like LSP 
with mass slightly below M\y can be a viable cold Dark Matter candidate, both from the 
phenomenological and from the model building point of view. 
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Appendix: Expressions for C —Functions 


In Sec. 2 we gave general expressions for the one-loop corrections from Yukawa interactions to 
the ZxiXi vertex, eqs.(12), as well as for the (pXiXi and A\?Xi couplings, eqs.(16), in terms of 
the C —functions defined in ref.[[L6[]. However, as already mentioned in Sec. 2, the expressions 
for the C —functions contain apparent divergencies both in the limit s —> 4m|o relevant for 
the calculation of the LSP relic density, and in the limit s —> 0 relevant for the calculation of 
the LSP-nucleon scattering cross section.]^] In case of the higher C —functions the necessary 
cancellations become too delicate for a reliable numerical treatment even if “double precision 1 ' 
variables are used. We have therefore re-evaluated the relevant Feynman parameter integrals 
in these two kinematical limits, which allows us to express the C—functions appearing in Sec. 2 
as combinations of B— functions. These expressions are collected in this Appendix. 

In our notation the scalar three-point function Cq is defined as 


G 0 (s, m 2 , Mi, M 2 , M 3 ) = — J dy jdx m 2 y 2 + s(x 2 — xy ) + y (m% — M, 


<1 


m 


^Recall that eqs.(12) and (16) have been written in a convention where both momenta k\ and point 
towards the vertex. In case of LSP-nucleon scattering the sign of one these momenta therefore has to be 
inverted. 
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+x(M 1 2 - Ml) + Ml - ie 


i -1 


(A.l) 


this definition coincides with that used in Appendix C of ref. [ 161. This gives: 


C 0 (4m 2 , m 2 , M 1? M 2 , M 3 ) = 1 \B 0 {m 2 , M 1; M 3 ) + B 0 (m 2 , M 2 , M 3 ) 


£> L 


C*o(0, m 2 , M l5 M 2 , M 3 ) — 


-2H 0 (4m 2 , M l5 M 2 
1 


(A. 2 a) 


Ml - Mf L 


Bo{in , M l5 M 3 ) — -B 0 (m 2 , M 2 , M 3 ) , (A.2b) 


where 


D = 2 (m 2 + M 3 2 ) - Ml - M 2 . 
In the limit —> M 2 , eq.(|A.2b|) reduces to: 


C 0 (0,m 2 ,M,M,M 3 ) = 


2 m 2 


, M 2 Ml + m 2 -M 2 

log —j + —-—- L 


Ml 


(A.3) 


(A.4) 


where we have introduced 

A = 2m 2 (M 2 + Ml) - m 4 - (ihf 2 - M 3 2 ) 2 ; 


L = 


2 arctan Af2+ )g_ m2 , 

, _ M 2 +M^—?n 2 + v / --A 
lQ g M 2 +M|—in 2 —\/—A ’ 


A > 0 
A < 0 


(A.5a) 
(A.5b) 


If eqs. (|A.2|) are used for C 0 , the function Cy defined in ref. [[H|] has apparent divergencies 
only at s —» 4m 2 (we suppress the imaginary infinitesimal — ie from now on): 


Ci (Am 2 , m 2 , Mi, M 2 , M 3 ) = — [ dy [ dx 

Jo Jo 


y/2 


o m 2 (y 2x ) 2 + y (Ml — Ml — m 2 ) + x (Mf — Mf) + Mf 

1 1 
AD + D 2 


{M| [S 0 (m 2 , Mi, M 3 ) + B 0 (m 2 , M 2 , M 3 ) 


+m 


B 2 (m 2 , Ms, Mi) + B 2 (m 2 , Ms, M 2 ) 


M 2 — /If 2 r 

+ 1 2 2 [Bi(m 2 , Ms, Mi) - Bi(m 2 , M 3 , M 2 ) 

- 2 (m 2 + Ml) + Mf - Mf] 5 0 (4m 2 , Mi, M 2 ) (A. 6 ) 

+2 (Mf - Mf) 5!(4m 2 , M 2 , Mi) + 8m 2 H 3 (4m 2 , Mi, M 2 )} , 

where D has been defined in eq.( |A.3| )- Here we have used the higher B functions as defined 
in ref.j44|; recall that our definition of Bi differs by an overall sign from that of ref. ||T 6 | . 


Similarly, after application of eqs. (|A.2| ), the function Ci contains apparent divergencies 
only at s = 0 : 


Ci (0, m 2 , Mi, M 2 , Ms) = — f dy [ dx 

Jo Jo 


x - y/2 


o y Jo m 2 y 2 + y (Mf - Ml - m 2 ) + x (M 2 - Ml) + Ml 
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2 {Ml - Ml 


1 + 2 


.Ml-Ml 

Ml - Mi 


B\ (m 2 , M 3 , Ml) - Bi (m 2 , M 3 , Ml) 


2 m 2 


Ml - Mi L 


B 3 (m 2 , M 3 , M 2 ) - B 3 {m 2 , M 3 , Ml) (A.7) 


2 Ml 


Ml - Mi L 


Bo{m 2 , M 2 , M 3 ) - B 0 {m 2 , Mi, M 3 ) 


- 1 


Note that by construction |l6| . Cj —> 0 as M 3 —> M 2 . 

The functions and Cl only appear in the Z vertex, eqs.(12). Moreover, the coefficient 
in front of Cl vanishes for s —> 4m|o, while the coefficient in front of Cl vanishes for s — > 0. 


We therefore only need to consider C 2 in the limit s —> 0: 

C^(0, m 2 , Mi, M 2 , M 3 ) — — f dy f dx 

Jo Jo 


2/ 2 /4 


o y Jo mV + y {Mi - Mf - m 2 ) + x {Ml - Mi) + Mi 
1 


5 2 (m , M 3 , M 2 ) - B 2 {m , M 3 , Mi) 


(A.8) 


4 {Ml - Mi) L 
In the limit M\ —> M 2 , this reduces to 

(7 2 + (0, m 2 , M, M, M 3 ) = 1 [^(m 2 , M 3 , M) + (m 2 , M 3 , M) + C 0 (0, m 2 , M, M, M 3 )] , (A.9) 

where C 0 (0, m 2 , M, M, M 3 ) is given in eq.( |A.4|) , and B' 0 and B[ are the derivatives of B 0 and 
Bi with respect to their first argunrent|j] Similarly, we need C 2 only in the limit s —> 4m 2 : 


C 2 (4m , m , Mi, M 2 , M 3 ) — — dy dx 


(x - y/2f 


1 

4D L 


0 u Jo m 2 {y — 2x) 2 + y {Mi — Mi — m 2 ) + x (Mf — Mi) + Mi 
B 2 {m 2 , M 3 , Mi) + H 2 (m 2 , M 3 , M 2 ) 

+8-B 3 (4m 2 , Mi, M 2 ) — 25 0 (4m 2 , Mi, M 2 ) , (A.10) 


where 14 is again given by eq.( |A-3|) . 

Finally, we note that the divergent function that appears in eqs.(12) can be computed 
from the general expression given in eq.(C4) of ref.JT6[, using the results for C 0 , Cf and C'f 
collected in this Appendix. 
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Figure Captions 


Fig. 2 : The chargino-LSP mass difference (solid), the axial-vector ZxiXi coupling (long 
dashed), and the h°XiXi coupling (short dashed) as a function of the soft breaking A param¬ 
eter, including one-loop corrections involving Yukawa couplings. “Low” and “high” leading- 
order estimates for the branching ratio for inclusive radiative b decays are also shown; for 
our assumption of exactly universal weak-scale soft breaking squark masses, this gives the 
strongest constraints on A in the region of interest, as discussed in the text. All quantities 
have been re-scaled, as indicated. 

Fig. 3 : The chargino-LSP mass difference (a), the LSP relic density Q^h 2 (b), and the 
expected LSP detection rate in a 76 Ge detector (c), as a function of the gaugino fraction, 
defined as the sum of the squares of the gaugino components of the LSP eigenvector. These 
results are for a fixed LSP mass, so that both M 2 and p, vary along the text. Further, |A| has 
been decreased from 2.7 rriq to 2.5rriq as M 2 was increased from about 150 GeV to 1 TeV. 

Fig. 4 : As in Fig. 3, but for positive sign of /q also, |A| has been kept fixed in this figure, 
and m^o has been reduced by 2 GeV, in order to be closer to the region of parameter space 
where the prediction for Q^h 2 is maximized. 

Fig. 5: In the region of the (q < 0, M 2 ) half-plane corresponding to a higgsino-like LSP, we 
show contours of constant Q^h 2 (dashed) and contours of constant LSP detection rate(CR) in 
a ' 6 Ge detector, in units of events/(Kg • day) (dotted and dot-dashed); note that the values 
on the latter contours are different for Figs, (a) (no Yukawa loop corrections), (b) (A < 0) 
and (c) (A > 0). The region to the right of the solid line is excluded by our interpretation of 
the LEP chargino search limit, as discussed in the text. 
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